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ABSTRACT 


Brown,  M.  W.  1985.  A  Fortran  Program  to  Calculate  Several 
Diversity  Indices.  U.S.  Department  of  Agriculture, 
Agricultural  Research  Service,  ARS-27,  2h  p. 

An  interactive  fortran  program  that  will  calculate 
Berger-Parker  dominance  index,  species  richness,  Simpson 
index,  Probability  of  Interspecific  Encounter,  Shannon-Weaver 
index,  and  Williams  index  is  presented.  The  program  also 
calculates  the  fit  to  the  log-series  distribution.  The 
concept  of  diversity  and  the  uses  and  limitations  of  each 
index  are  discussed.  A  sample  run  of  the  program  is  included 
with  comments  on  the  input  and  output. 

Keywords:  diversity,  Berger-Parker  dominance  index,  species 
richness,  Simpson  index,  Probability  of  Interspecific 
Encounter,  Shannon-Weaver  index,  Williams  index,  log-series 


Copies  of  this  publication  can  be  purchased  from  the  National 
Technical  Information  Service,  5285  Port  Royal  Road, 
Springfield,  Virginia  22161. 
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A  FORTRAN  PROGRAM  TO  CALCULATE  SEVERAL  DIVERSITY  INDICES1 
By  M.  W.  Brown2 


Diversity  is  a  property  of  a  collection  of  individuals  that 
can  be  classified  into  mutually  exclusive  groups.  Much  of  the 
theory  and  many  indices  were  initially  developed  for  use  in 
information  theory  (Yule  1944;  Shannon  and  Weaver  1949). 
Ecologists  have  adapted  this  information  theory  concept  of 
diversity  to  investigate  and  describe  the  structure  and 
function  of  communities  in  which  individuals  can  be  classified 
into  species  or  other  taxonomic  and/or  ecological  groupings. 
Some  of  the  same  indices  developed  for  information  theory  have 
been  used  in  community  ecology,  but  many  other  indices  have 
been  developed  as  well  and  the  theoretical  basis  has  been 
broadened  (Fisher  et  al.  1943;  Kempton  and  Taylor  1974;  May 
1975;  Pielou  1977). 

Diversity  is  generally  considered  to  be  composed  of  two 
elements:  number  of  species  or  categories,  and  the  relative 
abundances  of  those  species  (evenness).  Herein  lies  the 
difficulty  of  developing  a  universal  index  of  diversity;  that 
is,  incorporating  both  aspects  into  one  number.  Poole  (1974) 
stated  that  the  only  objective  measure  of  diversity  is  the 
number  of  species  (species  richness);  once  evenness  is  added, 
the  index  becomes  dependent  upon  how  evenness  is  defined. 
Another  hazard  in  combining  evenness  and  the  number  of  species 
in  a  single  index  is  that  their  individual  effects  may  be 
confounded  (King  and  Pimm  1983). 

Despite  these  difficulties,  diversity  indices  have  a  valuable 
role  in  community  ecology.  Because  of  different  assumptions 
and  definitions  of  evenness,  each  index  provides  a  slightly 
different  perspective  of  diversity.  One  must  keep  in  mind  the 
drawbacks  of  each  index  and  select  the  index  that  best  suits 
the  task  for  which  it  is  to  be  used.  The  researcher  must 
first  decide  what  he  or  she  wants  to  know  about  the  community 
and  why  (for  example,  to  compare  with  other  communities,  to 
compare  the  same  community  over  time,  to  evaluate  the  effect 
of  a  perturbation,  or  to  analyze  the  structure  or  functioning 
of  a  particular  community),  and  then  select  the  most 
appropriate  index.  Southwood  (1978)  proposed  a  possible 
approach  to  selecting  the  proper  index.  For  a  general  review 
of  diversity  and  several  of  the  indices,  see  Peet  (1974), 

Poole  (1974),  Pielou  (1977),  or  Southwood  (1978). 


1This  material  is  based  upon  work  supported  by  the  U.S. 
Department  of  Agriculture  under  Agreement  No.  58-32U4-4-675. 


2Adjunct  Assistant  Professor,  West  Virginia  University, 
Appalachian  Fruit  Research  Station,  Route  2,  Box  45, 
Kearneysville,  WV  25430 
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SPECIES  RICHNESS 


BERGER-PARKER 
DOMINANCE  INDEX 
(BERGER  and 
PARKER  1970) 


The  program  presented  here  calculates  a  number  of 
nonparametric  diversity  indices  and  one  parametric  index, 
Williams  index.  The  nonparametric  indices  are  not  dependent 
upon  the  underlying  distribution  of  the  data,  whereas  Williams 
index  is  based  on  the  log-series  distribution.  Other  proposed 
distributions,  such  as  the  geometric  (niche  preemption  model) 
(Poole  1974),  broken  stick  model  (MacArthur  1957),  and 
log-normal  (Preston  1948),  are  not  used  because  the  log-series 
seems  to  be  more  widely  applicable  and  has  a  broader 
theoretical  base  (Taylor  et  al.  1976).  The  nonparametric 
indices  calculated  are  species  richness,  Berger-Parker 
dominance  index,  Simpson  index,  Probability  of  Interspecific 
Encounter  and  Shannon-Weaver  index. 

The  eguations  used  assume  the  data  are  from  a  completely 
random  sample  of  the  community,  containing  all  species  which 
occur  in  the  community.  These  are  guite  rigid  assumptions, 
which  unfortunately  are  seldom  met.  Although  violation  of  the 
assumptions  will  bias  the  indices,  they  can  still  be  useful  as 
relative  estimates  of  diversity  if  the  limitations  are  kept  in 
mind. 

First,  I  present  a  discussion  about  each  index,  eguations 
used,  and  references  to  each.  The  Fortran  listing  (from  a 
Wang  VS)  is  then  given  with  an  example  and  explanations  of  the 
input  and  output. 

This  index  is  the  number  of  species  (S)  in  the  sample.  It  is 
only  one  part  of  what  constitutes  diversity,  omitting  evenness 
of  abundances.  However,  it  is  a  completely  objective  index  as 
long  as  the  sample  is  random  and  contains  all  species  in  the 
community . 

This  index  is  only  a  measure  of  evenness  with  no  information 
on  the  number  of  species  present.  The  Berger-Parker  index 
is  defined  as  the  probability  that  an  individual  selected  at 
random  from  the  community  belongs  to  the  most  abundant  species 
or,  eauivalently ,  the  proportion  of  all  individuals  that  are 
of  the  most  abundant  species.  It  is  calculated  as  the  largest 
Pi  where  Pi=nj/N,  N=the  total  number  of  individuals  in 
the  sample  and  n±=the  number  of  individuals  in  the  sample 
belonging  to  species  i.  The  index,  which  is  easy  to  calculate 
and  objective,  is  free  of  sample  bias  (Taylor  et  al.  1976), 
and  characterizes  the  distribution  of  species  abundances  as 
well  or  better  than  other  indices  (May  1975). 
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SIMPSON  INDEX 
(SIMPSON  1949) 


The  Simpson  index  (C)  is  the  probability  that  if  two 
individuals  are  drawn  at  random  from  a  community  they  will  be 
of  the  same  species.  It  is  calculated  as: 


E [n. (n. -1 )/N(N-l ) ] 


and  has  a  variance  of: 

4N(N-1) (N-2)Zp.  3  +  2N ( N— 1 ) Ep . 2  -  2N ( N-l ) ( 2N-3 ) ( Ep. 2 ) 2 


[ N ( N— 1 ) ]2 


This  index,  like  the  Berger-Parker,  is  a  measure  of 
concentration  or  dominance.  It  is  most  influenced  by  the 
relative  degree  of  abundance  of  the  few  most  abundant 
(dominant)  species  rather  than  by  the  overall  evenness  of  all 
species.  Because  it  measures  dominance,  as  diversity 
increases,  C  decreases.  To  transform  this  to  an  index  (D) 
that  is  proportional  to  diversity,  three  conversion  eguations 
have  been  proposed.  Pielou  (1977)  suggested  D=  -In  (C)  in 
which  form  it  is  comparable  to  the  Shannon-Weaver  index  and  is 
a  special  case  of  a  general  diversity  function  given  by  Renyi 
(1961).  Berger  and  Parker  (1970)  and  Hurlbert  (1971) 
suggested  D=l-C  where  D  becomes  the  probability  of  the  next 
encounter  (the  next  individual  sampled)  being  a  different 
species.  Finally,  Southwood  (1978)  suggested  D=l/C.  All 
three  forms  of  D  are  calculated  by  the  program. 

Species  evenness,  defined  as  the  ratio  of  observed  diversity 
to  maximum  diversity  (that  is,  all  species  in  the  sample  being 
egually  abundant),  can  be  a  useful  measure.  Evenness  can 
range  from  a  value  of  0.0  if  the  observed  diversity  is  at  a 
minimum  (that  is,  one  species  with  all  the  individuals  except 
for  one  individual  in  each  of  the  other  species)  to  1.0  if  the 
observed  diversity  is  the  maximum  possible  for  the  sample. 
Through  the  calculation  of  evenness,  the  diversity  index  is 
standardized,  allowing  for  more  valid  comparisons  among 
communities  or  among  samples  from  the  same  community. 

Evenness,  based  on  the  Simpson  index  is  calculated  as: 


c 


(s-i)/s 


max 


C  .  =  (2N-S) (S-l)/N 

min 


Evenness  =  (C 


max 
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PROBABILITY  of 
INTERSPECIFIC 
ENCOUNTER  (PIE) 
(HURLBERT  1971) 


SHANNON-WEAVER 
INDEX  (SHANNON 
AND  WEAVER  1949) 


This  index  is  very  similar  to  the  Simpson  index.  PIE  is 
defined  as  the  proportion  of  all  the  possible  encounters 
among  individuals  that  are  interspecific.  The  difference 
between  the  Simpson  index  and  PIE  is  the  same  as  sampling  from 
a  finite  population  with  replacement  (Simpson  index)  or 
without  replacement  (PIE).  For  all  but  very  small  data  sets 
PIE  =  1-Simpson  index,  but  they  are  different  conceptually. 

PIE  is  calculated  as: 


£ ( n . /N ) [ (N-n. )/ (N— 1 ) ] 


The  variance  is  calculated  as  the  variance  of  the  Simpson 
index  multiplied  by  [N/(N-1)]2.  Species  evenness  can  be 
defined  as  (PIE  -  PIEmj_n)/(PIE  max  -  PIE  mj_n)  with  the 
following  formulas: 

pie  =  [N/(N-1)][(S-1)/S] 
max 

PIE  .  =  [N/(N— !)][( 2N— S ) ( S— 1 ) /N  2 ] 


This  index  is  the  most  commonly  used  and  accepted  in 
ecological  work.  It  has  been  criticized,  however,  as  having 
no  intuitive  ecological  meaning  (Taylor  et  al.  1976). 

Southwood  (1978)  considered  it  to  be  a  distraction  rather  than 
an  asset  in  ecological  use.  It  is  easily  calculated,  and  it 
is  recognized  by  most,  if  not  all,  ecologists.  The 
Shannon-Weaver  index  (H')  is  a  measure  of  uncertainty  in 
predicting  the  species  of  an  individual  drawn  at  random  from 
the  community.  Two  eguations  are  used  by  the  program.  The 
first : 


-Zp^ [ln(pd  ] 

is  the  one  most  used  but  is  biased.  The  second  is  the 
expected  value  of  H'.  Only  the  first  four  terms  of  the 
expansion  given  by  Hutcheson  (1970)  are  used  in  the  program. 


-Zp.[ln(p.)]  -  (  S— 1 )  /  2N  +  (1-Zp.'1  )/12N2  +  Z  (  pD1 -p  “2  ) /12N 3 

The  variance  of  E(H'),  using  the  first  three  terms  of 
Hutcheson's  (1970)  expansion,  is: 

(Zp.[ln(pd]2  -  {Zp.  [  ln(p.)  ]  } 2  )/N  +  (S-1)/2N2  + 

(-1  +  Zp  F1  -  {  Zp/1  [  ln(pd  ]  }  +  Zp.'1  {Zpd  ln(pd  ]  }  )/6N3 
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WILLIAMS  INDEX 
(FISHER  ET  AL. 
1943) 


The  maximum  diversity  for  a  given  community  is  H'max=ln  (S). 
This  can  be  used  to  estimate  evenness  of  the  community, 
H'/H'max*  Strictly  speaking,  the  index  and  its  variance  are 
inestimable  when  the  true  number  of  species  in  the  community 
is  unknown. 

Williams  index  is  the  only  parametric  diversity  index 
calculated  by  this  program,  the  index  being  a  parameter  (  a  ) 
of  the  log-series  distribution.  A  justification  for  the  use 
of  this  index  is  given  by  Kempton  and  Taylor  (1974)  and  Taylor 
et  al.  (1976).  They  stated  that  this  index  should  be  the 
index  of  choice,  assuming  a  log-series  distribution  a  priori. 
The  log-series,  an  approximation  to  the  Gamma  distribution, 
has  two  parameters:  a  (Williams  index),  and  X.  Alpha  is  a 
characteristic  of  the  community.  The  value  of  a  is  most 
influenced  by  the  moderately  abundant  species.  Therefore,  it 
is  relatively  insensitive  to  rare  or  dominant  species,  the 
presence  and  abundance  of  which  are  more  subject  to  sampling 
errors.  X  is  a  characteristic  of  the  sample  itself  and  is 
dependent  on  the  total  number  of  individuals. 

Alpha  is  calculated  by  iteration  of  the  maximum  likelihood 
estimator  (Taylor  et  al.  1976): 

S  =  a [ In ( 1+N/a )  ] 

X  is  calculated  by  solving  eguations  (Pielou  1977): 

S  =  -a [ In ( 1-X ) ] 


and 


N  =  aX/ ( 1— X ) 


fora,  and  iterating  the  result: 

S  =  { [ N ( l-X ) ] /X }  [ -In ( 1— X ) ] 

The  program  prints  out  the  predicted  number  of  species  from 
each  iteration.  These  values  can  be  used  to  check  for  errors 
if  the  predicted  S  is  very  different  from  the  actual  S. 
Another  check  can  be  made  by  using  figure  8  of  Fisher  et  al. 
(1943).  Variance  of  a  is  calculated  by 

-a/ln ( 1-X ) 
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Calculation  of  the  expected  log-series  distribution  using  the 
calculated  parameters  a  and  X  is  by 

E ( r )  =  a xr/ r 

where  E(r)  is  the  expected  number  of  species  with  r 
individuals.  A  chi-square  goodness  of  fit  test  is  conducted 
on  the  data,  combined  into  octaves  (see  Poole  1974),  to  judge 
whether  a  should  be  used  as  the  diversity  index.  Finally, 
expected  values  of  S  and  N  based  on  the  log-series 
distribution  are  calculated  as: 

S  =  -a [ In ( 1— X ) ] 

N  =  aX/(l-X) 
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APPENDIX  A:  PROGRAM  INPUT 


Note:  Computer  generated  comments  and  questions  are  in 
capital  letters;  user  responses  are  in  lowercased  letters  or 
numerals.  Comments  by  the  author  are  offset  with  an  asterisk. 

INPUT  TITLE  (UP  TO  32  CHARACTERS) 

abandoned  orchard  wva,  18-iv-83 

*Data  on  relative  abundance  of  herbivorous  insects  found  in  an 
abandoned  apple  orchard  in  West  Virginia  collected  by  M.  W. 
Brown  and  Cynthia  R.  Loerch,  Appalachian  Fruit  Research 
Station,  Kearneysville,  WV. 

INPUT  NUMBER  OF  QUADRATS.  (TWO  DIGITS) 

06 

*Zeros  are  put  in  to  show  spacing  for  the  proper  format  and 
are  not  needed  when  supplying  data. 

INPUT  TOTAL  NUMBER  OF  SPECIES.  (TWO  DIGITS) 

12 

*This  is  the  number  of  species  found  on  all  quadrats  in  the 
sample. 

BEGIN  DATA  ENTRY. 

INPUT  DATA  FOR  SPECIES  1 
000700010002000000040003 

*Zeros  are  included  only  to  show  spacing.  There  are  6  fields 
(number  of  fields  taken  from  the  number  of  quadrats  supplied 
earlier)  of  4  digits  (F4.0  format). 

INPUT  DATA  FOR  SPECIES  2 

000300030006000600080007 

INPUT  DATA  FOR  SPECIES  3 

000400050005000700020000 

INPUT  DATA  FOR  SPECIES  4 

001 000000001000000000000 


ll 


INPUT  DATA  FOR  SPECIES  5 


OOOOOOOOOOOIOOOCOOOOOOOO 
INPUT  DATA  FOR  SPECIES  6 
000000000003000200030004 
INPUT  DATA  FOR  SPECIES  7 
0000000000000001 0001 0000 
INPUT  DATA  FOR  SPECIES  8 
000000000000000100000000 
INPUT  DATA  FOR  SPECIES  9 
0000000000000001 00000000 
INPUT  DATA  FOR  SPECIES  10 
000000000000000000030005 
INPUT  DATA  FOR  SPECIES  11 
000000000000000000000001 
INPUT  DATA  FOR  SPECIES  12 
000000000000000000000001 
DO  YOU  WANT  TO  CHECK  DATA?  Y/N 

y 

SPECIES  QUADRAT 


1 

2 

3 

4 

5 

6 

1 

7. 

1. 

2. 

0. 

4. 

3 

2 

3. 

3. 

6. 

6. 

8. 

7 

3 

4. 

5. 

5. 

7. 

2. 

0 

4 

o’ 

• — 1 

0. 

1. 

0. 

0. 

0 

5 

0. 

0. 

1. 

0. 

0. 

0 

12 


6 

0. 

0. 

3. 

2. 

3. 

4. 

7 

0. 

0. 

0. 

1. 

1. 

0. 

8 

0. 

0. 

0. 

1. 

0. 

0. 

9 

0. 

0. 

0. 

1. 

0. 

0. 

10 

0. 

0. 

0. 

0. 

3. 

5. 

11 

0. 

0. 

0. 

0. 

0. 

1. 

12 

0. 

0. 

0. 

0. 

0. 

1. 

DO  YOU  WANT  TO  CORRECT  ANY  DATA?  Y/N 

y 

WHICH  QUADRAT?  (TWO  DIGITS) 

01 

WHICH  SPECIES?  (TWO  DIGITS) 

04 

ENTER  DATA  FOR  QUADRAT  1  SPECIES  4 
0000 

DO  YOU  WANT  TO  CORRECT  ANY  OTHER  DATA?  Y/N 
n 

DO  YOU  WANT  TO  LOOK  AT  THE  DATA  AGAIN?  Y/N 
n 

*At  this  point  you  could  examine  the  data  again  and  correct 
more  errors  by  answering  "Y". 

DIVERSITY  CALCULATIONS  IN  PROGRESS 
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APPENDIX  B:  PROGRAM  OUTRJT 


ABANDONED  ORCHARD  WVA,  18-IV-83 


********** 

DATA 

********** 

SPECIES 

QUADRAT 

1 

2  3 

4 

3  6 

1 

7. 

1. 

2. 

0. 

4. 

3. 

2 

3. 

• 

6. 

6. 

8. 

7. 

3 

4. 

5. 

5. 

7. 

2. 

0. 

4 

0. 

0. 

1. 

0. 

0. 

0. 

3 

0. 

0. 

1. 

0. 

0. 

0. 

6 

0. 

0. 

3. 

2. 

3. 

4. 

7 

0. 

0. 

0. 

1. 

1. 

0. 

8 

0. 

0. 

0. 

1. 

0. 

0. 

9 

0. 

0. 

0. 

1. 

0. 

0. 

10 

0. 

0. 

0. 

0. 

3. 

5. 

11 

0. 

0. 

0. 

0. 

0. 

1. 

12 

0. 

0. 

0. 

0. 

0. 

1. 

Ih 


*****  DIVERSITY  statistics  ***** 


SPECIES 

NUMBER 

PROBABILITY 

1 

17. 

0.1683 

2 

33. 

0.3267 

3 

23. 

0.2277 

4 

1. 

0.0099 

5 

1. 

0.0099 

6 

12. 

0.1188 

7 

2. 

0.0198 

8 

1. 

0.0099 

9 

1. 

0.0099 

10 

8. 

0.0792 

11 

1. 

0.0099 

12 

1. 

0.0099 

SPECIES 

RICHNESS  =  12 

101.  INDIVIDUALS 

BERGER-PARKER  DOMINANCE  INDEX 

=  0.3267 

SIMPSON  INDEX  OF  DIVERSITY  =  0.2004  VARIANCE  =  0.000435 


1 /SIMPSON  INDEX  = 

4.9901 

1-SIMPSON  INDEX  = 

0.7996 

-LN  SIMPSON  INDEX  = 

1.6075 

MINIMUM  SIMPSON  = 

0.0742 

MAXIMUM  SIMPSON  = 

0.7931 

EVENNESS  = 

0.8244 

PROBABILITY  OF  INTERSPECIFIC  ENCOUNTER  =  0.7996  VARIANCE  =  0.000443 


MAXIMUM  P.I.E.  =  0.9258 
MINIMUM  P.I.E.  =  0.2069 
EVENNESS  =  0.8244 


SHANNON-WEAVER  INDEX  H=  1.8081 

EXPECTED  SHANNON-WEAVER  INDEX  =  1.7429  VARIANCE  =  0.008775 

MAXIMUM  DIVERSITY  =  2.4849 

EVENNESS  =  0.7276 
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*****  LOG-SERIES  STATISTICS  ***** 


ALPHA  (WILLIAMS  INDEX)  =  3.55  (PREDICTED  NUMBER  OF  SPECIES  = 

X  =  0.9661  (PREDICTED  NUMBER  OF  SPECIES  = 

VARIANCE  OF  ALPHA  =  1.048949 


CHI-SQUARE  GOODNESS  OF  FIT  TABLE 


INDIVIDUALS  PER 
(OCTAVE) 

SPECIES  NUMBER 

EXPECTED 

OF  SPECIES 

OBSERVED 

CHI-SQUARE 

1 

3.430 

6 

1.9263 

2-3 

2.724 

1 

1 . 0909 

4-7 

2.250 

0 

2.2501 

8-15 

1.766 

2 

0.0311 

16-31 

1.170 

2 

0.5885 

32-63 

0.541 

1 

0.3890 

64+ 

0.125 

0 

0.1255 

GOODNESS  OF  FIT 

TO  THE  LOG-SERIES  MODEL 

CHI-SQUARE  = 

6.4014 

EXPECTED  NUMBER  OF  SPECIES  =  12.0144 

EXPECTED  NUMBER  OF  INDIVIDUALS  =  101.1695 


12.0086  ) 
11.9943  ) 
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APPENDIX  C:  PROGRAM  LISTING 


C  JOB  CONTROL  STATEMENTS  MUST  BE  PLACED  HERE 
C 

C  PROGRAM  TO  READ  DATA  ON  SPECIES  ABUNDANCE  PER  QUADRAT  AND 
C  CALCULATE  DIVERSITY  INDICES  AND  FIT  TO  THE  LOG-SERIES  MODEL. 

C 

C  PROGRAM  VARIABLES  ARE  DIMENSIONED  TO  HANDLE  UP  TO  100  SPECIES  AND 
C  10  QUADRATS  BUT  IS  FORMATTED  FOR  10  SPECIES  WITH  AN  ABUNDANCE  OF 
C  UP  TO  9999.  INDIVIDUALS.  THESE  RESTRICTIONS  ARE  ONLY  FOR  FORMATTING 
C  AND  MAY  BE  ALTERED  BY  EDITING  THE  CORRESPONDING  DIMENSION  AND 
C  FORMAT  STATEMENTS. 

C 

INTEGER  Q( 10) ,OBS( 7) 

REAL  N(100,10) ,SPLIST(100) ,P(100) , EXPECT (7) ,CHI(7) 

REAL*8  A,B,C,D,E 

DATA  N/1000*0. 0/ , BERG/O. 0/ ,SPLIST/100*0. 0/ , TOTAL/O. 0/ ,RRRP/0. 0/ 
DATA  0BS/7*0/ , SIMP/O . 0/ , SHAM/O . 0/ , RP/O. 0/ , PSQ/O. 0/ , RRP/O . 0/ 

DATA  PIE/O . 0/ , SUMP3/0 . 0/ , SUMP2/0. 0/ , VARSIM/O . 0/ 

WRITE(0, 900) 

900  FORMAT ( '  INPUT  TITLE  (UP  TO  32  CHARACTERS)') 

READ (0, 902)  A,B,C,D 
902  FORMAT (4A8) 

WRITE(6,904)  A,B,C,D 
904  FORMAT ( ' 1 ' , 5X , 4A8 ) 

WRITE( 0, 906) 

906  FORMAT ( '  INPUT  NUMBER  OF  QUADRATS.  (TWO  DIGITS)') 

READ(0, 908)  J 
908  FORMAT (12) 

WRITE(0,910) 

910  FORMAT ( '  INPUT  TOTAL  NUMBER  OF  SPECIES.  (TWO  DIGITS)') 

READ(0, 912)  I 
912  FORMAT ( 12) 

S=I 

WRITE(0, 914) 

914  FORMAT ( '  BEGIN  DATA  ENTRY.') 

DO  100  JJ=1,J 
Q( JJ)=JJ 
100  CONTINUE 

DO  102  11=1,1 
WRITE(0,916)  II 

916  FORMAT ( '  INPUT  DATA  FOR  SPECIES', 13) 

READ( 0, 918)  (N( II , JJ) , JJ=1 , J) 

918  FORMAT ( 10F4. 0) 

102  CONTINUE 

WRITE (0,920) 

920  FORMAT ( '  DO  YOU  WANT  TO  CHECK  DATA?  Y/N') 

READ( 0,922)  E 
922  FORMAT (A8) 

IF(E.NE.'Y')  GO  TO  304 
300  WRITE(0,924) 

924  FORMAT(////, '  SPECIES' ,10X, 'QUADRAT' ) 

WRITE(0, 926)  (Q( JJ) , JJ=1, J) 

926  F0RMAT(//,9X,10(I2,3X) ,//) 


IT 


DO  104  11=1,1 

WRITE(0, 928)  II , (N( II , JJ) , JJ=1 , J) 

928  FORMAT ( 15, 2X, 10F5. 0) 

IF  (II.EQ.20)  PAUSE 
IF  (II.EQ.40)  PAUSE 
IF  (II.EQ.60)  PAUSE 
104  CONTINUE 

WRITE(0,930) 

930  FORMAT ('  DO  YOU  WANT  TO  CORRECT  ANY  DATA?  Y/N') 

READ(0, 922)  E 
IF(E.NE. ' Y ' )  GO  TO  304 
302  WRITE(0, 932) 

932  FORMAT ( '  WHICH  QUADRAT?  (TWO  DIGITS)') 

READ (0,934)  JJ 
934  FORMAT ( 12) 

WRITE( 0, 936) 

936  FORMAT ( '  WHICH  SPECIES?  (TWO  DIGITS)') 

READ( 0, 938)  II 
938  FORMAT ( 12) 

WRITE(0,940)  JJ, II 

940  FORMAT ( '  ENTER  DATA  FOR  QUADRAT* , 12, '  SPECIES', 13) 

READ(0, 942)  N(II,JJ) 

942  FORMAT ( F4. 0) 

WRITE( 0, 944) 

944  FORMAT ('  DO  YOU  WANT  TO  CORRECT  ANY  OTHER  DATA?  Y/N') 

RE AD (0,922)  E 
IF(E.EQ.'Y')  GO  TO  302 
WRITE(0,946) 

946  FORMAT ( '  DO  YOU  WANT  TO  LOOK  AT  THE  DATA  AGAIN?  Y/N') 

READ(0, 922)  E 
IF(E.EQ.’Y')  GO  TO  300 
304  WRITE(6,948) 

948  FORMAT( /////, 25X, ***********  DATA  ********** « f /// f •  SPECIES', 
120X, 'QUADRAT' ) 

WRITE(6, 950)  (Q(JJJ), JJJ=1,J) 

950  F0RMAT(//,10X,18I7) 

WRITE(6,952) 

952  FORMAT (/) 

DO  106  111=1,1 

WRITE(6, 954)  III , (N( III , JJJ) , JJJ=1 , J) 

954  FORMAT ( 15, 6X , 18F7. 0) 

106  CONTINUE 

WRITE(0, 956) 

956  FORMAT ( '1 ' , '  DIVERSITY  CALCULATIONS  IN  PROGRESS') 

DO  108  K=1 , I 
DO  110  KK=1,J 

SPLIST(K)=SPLIST(K)+N(K,KK) 

110  CONTINUE 
108  CONTINUE 

DO  112  K=1 , I 
TOTAL=TOTAL+SPLIST(K) 

112  CONTINUE 

DO  114  K=1 , I 
P(K)=SPLIST (K) /TOTAL 
SUMP2=SUMP2+P( K ) *P( K ) 

SUMP3=SUMP3+P(K)**3. 

RP=RP+(1./P(K)) 
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RRP=RRP+( (l./P(K))-(l./(P(K)*P(K)))) 

SIMP=SIMP+ ( SPL 1ST ( K )* ( SPLIST ( K ) -1 . ) ) / ( TOTAL*( TOTAL-1 . ) ) 
PIE=PIE+(SPLIST(K) /TOTAL )*( ( TOTAL-SPLIST(K) )/(TOTAL-l . ) ) 

IF(P(K) .GT.BERG)  BERG=P(K) 

IF(SPLIST(K) .LE.O.O)  GO  TO  114 
SHAN=SHAN-(P(K)*ALOG(P(K) ) ) 

PSQ=PSQ+( P ( K ) * ALOG ( P(K) ) *ALOG ( P ( K ) ) ) 

RRRP=RRRP+ ( (1 ./P(K) ) *ALOG (P(K) ) ) 

114  CONTINUE 

RSIMP1=1 . /SIMP 
RSIMP2=1 .-SIMP 
RSIMP3=-( ALOG( SIMP) ) 

VARSIM= ( ( 4 . *TOTAL* ( TOTAL-1 . ) * ( TOTAL-2 . ) *SUMP3 ) + ( 2 . *TOTAL* ( TOTAL-1 . 
1)*SUMP2)-(2.*T0TAL*(T0TAL-1 . )*( ( 2. *T0TAL)-3. )*(SUMP2*SUMP2) ) )/ 

2( ( TOTAL*(TOTAL-l . ) )*( TOTAL*( TOTAL-1 . ) ) ) 

VARPIE=( TOTAL/( TOTAL-1 . ) )*( TOTAL/( TOTAL-1 . ) )*VARSIM 
PIEMAX=( TOTAL /(TOTAL-1 . ) )*( (S-l . )/S) 

PIEMIN-( TOTAL/ (TOTAL-1 . ) )*( ( ( (2.*T0TAL)-S)*(S-1 . ) )/( TOTAL* TOTAL) ) 
PIEEVN=(PIE-PIEMIN)/(PIEMAX-PIEMIN) 

SIMMAX=( ( TOTAL-S+1 . )*(TOTAL-S) )/( TOTAL* (TOTAL-1 . ) ) 

SIMMIN=S*( ( (TOTAL/S)*( (TOTAL/S)-l . ) )/(TOTAL*( TOTAL-1 . ) ) ) 

SIMEVN=( SIMMAX-SIMP)/( SIMMAX-SIMMIN) 

SHANSQ=SHAN*SHAN 

ESHAN=SHAN-( ( 1-1 . )/(2*T0TAL) )+( (l.-RP)/(12.*T0TAL*T0TAL) )+ 
1(RRP/(12.*T0TAL**3. ) ) 

VARSHA= ( ( PSQ-SHANSQ ) /TOTAL ) + ( ( 1-1 . ) / ( 2 . *TOTAL*TOTAL ) ) + 

1 ( -1 . +RP-RRRP- ( RP*SHAN ) ) / ( 6 . *T0TAL**3 . ) 

HMAX=ALOG(S) 

EVEN=SHAN/HMAX 

WRITE(6,958) 

958  FORMAT ( '1* ,20X, ******  DIVERSITY  STATISTICS  *****',///, 

1'  SPECIES  NUMBER  PROBABILITY') 

DO  116  K=1 , I 

WRITE(6,960)  K,SPLIST(K) ,P(K) 

960  FORMAT ( '  ' , 17, 10X,F6.0,10X,F10.4) 

116  CONTINUE 

WRITE (6, 962)  I, TOTAL 

962  FORMAT ( ///, '  SPECIES  RICHNESS  =  ',I5,F10.0,'  INDIVIDUALS') 

WRITE (6, 964)  BERG 

964  FORMAT(///,'  BERGER-PARKER  DOMINANCE  INDEX  =  ',F10.4) 

WRITE( 6, 966)  SIMP,VARSIM,RSIMP1,RSIMP2,RSIMP3,SIMMIN,SIMMAX,SIMEVN 
966  FORMAT ( ///, '  SIMPSON  INDEX  OF  DIVERSITY  =  ' ,F15. 4, 5X, ' VARIANCE  =  ', 
1F10.6,//,6X, '  1/SIMPSON  INDEX  =  ' ,F9.4,/,6X, '  1-SIMPSON  INDEX  =  ', 
2F9.4,/,6X,'  -LN  SIMPSON  INDEX  =  ' ,F7.4,/,6X, '  MINIMUM  SIMPSON  = 
3F9.4,/,6X, '  MAXIMUM  SIMPSON  =  ' ,F9.4,/,6X, '  EVENNESS  =  *,F16.4) 
WRITE( 6 , 968)  PIE , VARPIE , PIEMAX , PIEMIN , PIEEVN 
968  FORMAT(///,'  PROBABILITY  OF  INTERSPECIFIC  ENCOUNTER  =',F15.4, 

15X, 'VARIANCE  =' ,F10.6,//,6X, '  MAXIMUM  P.I.E.  =' ,F12.4,/,6X, 

2'  MINIMUM  P.I.E.  =' ,F12.4,/,6X, '  EVENNESS  =',F18.4) 

WR I TE ( 6 , 970 )  SHAN , ESHAN , VARSHA , HMAX , E VEN 
970  FORMAT (///, '  SHANNON-WEAVER  INDEX  H=  ' ,F10.4,//,5X, '  EXPECTED 
1  SHANNON-WEAVER  INDEX  =' ,F10. 4, 5X ,' VARIANCE  =  ' ,F10.6,/,5X, 

2'  MAXIMUM  DIVERSITY  =  ' ,F9.4,/,5X, '  EVENNESS  =  ',F10.4) 

C  SECTION  TO  ESTIMATE  ALPHA  OF  THE  LOG-SERIES  DISTRIBUTION 
C  USES  MAXIMUM  LIKELIHOOD  ESTIMATOR 
C  S=ALPHA*(LOGE( 1+N/ALPHA) 

01=99999999. 


DO  118  K=l, lOOOO 
ALPHA=0 . l+K/100 . 

T=ALPHA*ALOG( 1 . +TOTAL/ALPHA ) 

D2=ABS(S-T) 

IF  (D2.GE.D1)  GO  TO  306 

Dl=D2 

XT=T 

XALPHA=ALPHA 
118  CONTINUE 
306  WRITE(6,972) 

972  FORMAT ( '  1 ' , 20X , ' *****  LOG-SERIES  STATISTICS  *****') 

WRITE(6, 974)  XALPHA ,XT 

974  FORMAT ( ///, '  ALPHA  (WILLIAMS  INDEX)  =  ' , Fll . 2 , '  (PREDICTED  NUMBER 
1  OF  SPECIES  =  ' ,F10.4, '  ) ' ) 

C  SECTION  TO  ESTIMATE  X  OF  THE  LOG-SERIES  DISTRIBUTION 
C  USES  MAXIMUM  LIKELIHOOD  ESTIMATOR 
C  S=(N( 1-X)/X)*(-L0GE( 1-X) ) 

T=0. 0 

Dl=99999999. 

X=0.1 

DO  120  K=1 , 8999 
X=0. l+K/10000. 

T=-( ( TOTAL* ( 1 .-X) )/X)*(ALOG( 1 .-X) ) 

D2=ABS(S-T) 

IF  (D2.GE.D1)  GO  TO  308 

Dl=D2 

XT=T 

XX=X 

120  CONTINUE 

308  WRITE(6,976)  XX, XT 

976  FORMAT ( /, '  X' ,21X, ,F12.4, '  (PREDICTED  NUMBER  OF  SPECIES  =' , 
1F10.4,'  )') 

VARA=-XALPHA/ALOG ( 1 . -XX ) 

WRITE(6, 978)  VARA 

978  FORMAT ( /  ,  1  VARIANCE  OF  ALPHA  =',F12.6) 

C  SECTION  TO  CALCULATE  FIT  TO  THE  LOG-SERIES  MODEL 
EXPECT ( 1 ) =XALPHA*XX 

EXPECT ( 2 ) = ( XALPHA* ( XX**2 . ) ) /2 . + ( XALPHA* ( XX**3 . ) ) /3 . 

EXPECT (3)=0. 

DO  122  K=1 , 4 

EXPECT ( 3)=EXPECT (3)+( XALPHA*( XX**(K+3. ) ) ) /(K+3. ) 

122  CONTINUE 

EXPECT (4)=0. 

DO  124  K=1 ,8 

EXPECT ( 4 ) =EXPECT ( 4 ) + ( XALPHA* ( XX** (K+7.)))/(k+7.) 

124  CONTINUE 

EXPECT ( 5)=0. 

DO  126  K=1 ,16 

EXPECT ( 5)=EXPECT ( 5)+( XALPHA*( XX**(K+15. ) ) )/(K+13. ) 

126  CONTINUE 

EXPECT (6)=0. 

DO  128  K=1 , 32 

EXPECT ( 6) =EXPECT (6)+(XALPHA*(XX**(K+31 . ) ) )/(K+31 . ) 

128  CONTINUE 

EXPECT (7)=0. 

DO  130  K=1 ,64 

EXPECT ( 7)=EXPECT ( 7)+( XALPHA* ( XX** (K+63. ) ) )/ (K+63. ) 
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o  o  o 


130  CONTINUE 

C  SECTION  TO  CALCULATE  FIT  TO  THE  LOG-SERIES  MODEL 
DO  132  M=1 , I 

IF  (SPLIST(M).EQ.l.)  0BS(1)=0BS(1)+1 
IF  (SPLIST (M) .EQ.2. .OR. SPLIST (M) .EQ.3. )  0bS(2)=0BS(2)+l 
IF  (SPLIST(M) . GE.4. . AND. SPLIST (M) .LT.8. )  0BS(3)^0BS(3)+1 
IF  ( SPLIST (M) .GE . 8. . AND. SPl 1ST (M) ,LT . 16. )  0BS(4)=0BS(4)+1 
IF  (SPLIST(M) .GE . 16. .AND . SPLIST(M) .LT . 32. )  OBS( 5)=0BS( 5) +1 
IF  ( SPLIST (M) .GE . 32. .AND. SPLIST (M) .LT . 64. )  0RS(6)=0BS(6)+1 
IF  (SPLIST(M) . GE.64. )  OBS( 7)=OBS(7)+l 
132  CONTINUE 
CHIS0=0 . 0 
WRITE (6, 980) 

980  FORMAT (////,20X, ' CHI -SQUARE  GOODNESS  OF  FIT  TABLE',///, 

1'  INDIVIDUALS  PER  SPECIES  NUMBER  OF  SPECIES  CHI- 

2SQUARE ' , /, 7X, ' (OCTAVE ) ' , 1 6X , 'EXPECTED  OBSERVED' , //) 

DO  134  M=l,7 

CHI (M)=(OBS(M) -EXPECT (M))*(OBS(M) -EXPECT (M) ) /EXPECT (M) 
CHISQ=CHISQ+CHI(M) 

134  CONTINUE 

WRITE(6,982)  (EXPECT(M) ,0BS(M),CHI(M),M=1,7) 

982  FORMAT ( '  ' , 8X, • 1 ' , 19X,F1C. 3, I10,F20.4, / , 8X, ' 2-3 ' , 18X,F10. 3, 110, 
1F20.4 , / , 8X, ' 4—7 ' , 18X,F10. 3, I10,F20. 4 , / , 8X, ' 8-15 ' , 17X, F10. 3 , 
2110, F20. 4, /, 7X, '16-31' ,17X,F10. 3,110, F20. 4,/, 7X, '32-63' , 
317X,F10. 3, I10,F20. 4, / , 7X, ' 64+ ' , 19X ,F10. 3, 110, F20. 4) 

WRITE( 6, 984)  CHISQ 

984  FORMAT( // , '  GOODNESS  OF  FIT  TO  THE  LOG-SERIES  MODEL', 

1 '  CHI-SQUARE  =' ,F12.4) 

EXPS=-XALPHA*ALOG( 1 . -XX ) 

EXPN=XALPHA*XX/ ( 1 . -XX ) 

WRITE(6,986)  EXPS,EXPN 

986  FORMAT ( // , '  EXPECTED  NUMBER  OF  SPECIES  =  ',F19.4,/, 

1'  EXPECTED  NUMBER  OF  INDIVIDUALS  =  ' ,F14.4) 

FOR  SOME  COMPILERS  A  "STOP"  STATEMENT  MAY  BE  NEEDED  HERE. 

END 
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FXt-HANGE  Rec’ST': 


